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It is shown that the generation linewidth of an auto-oscillator with a nonlinear frequency shift 
(i.e. an auto-oscillator in which frequency depends on the oscillation amplitude) is substantially 
larger than the linewidth of a conventional quasi-linear auto-oscillator due to the renormalization 
of the phase noise caused by the nonlinearity of the oscillation frequency. The developed theory, 
when applied to a spin-torque nano-contact auto-oscillator, predicts a minimum of the generation 
linewidth when the nano-contact is magnetized at a critical angle to its plane, corresponding to the 
minimum nonlinear frequency shift, in good agreement with recent experiments. 
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It is well-known that the linewidth Tq of a passive os- 
cillating circuit is determined by the ratio of its dissipa- 
tive element (e.g. resistance R) to its reactive element 
(e.g. inductance L): Tq = R/2L. When the oscillating 
circuit is connected to an active element (transistor, vac- 
uum tube, tunnel diode, etc.) and a source of a constant 
voltage (e.g. battery) the auto-generation of constant- 
amplitude oscillations at the resonance frequency of the 
oscillating circuit (uj = 1/%/LC, where C is the circuit 
capacitance) can take place [l], 0] • The equilibrium am- 
plitude of these auto-oscillations is determined by the 
dynamic balance between the positive nonlinear damping 
of the oscillating system and negative nonlinear damping 
introduced into the system by the active element [H, 0] • 

It is also well-established that the generation linewidth 
Au> in a typical auto-oscillator is determined, for the most 
part, by thermal phase noise (see e.g. Eq. (9.36) in [H), 
and can be expressed in the following general form, 



Auj = r 



E(a) ' 



(1) 



where fee is the Boltzmann constant, T is the absolute 
temperature, E(a) — (3\a\ 2 is the averaged energy of the 
auto-oscillation having the complex amplitude a, and (3 
is the coefficient relating the averaged energy to the auto- 
oscillation amplitude. For example, in an auto-oscillator 
with a standard linear oscillating circuit, [3 = C/2, where 
C is the capacitance of the oscillating circuit and a is the 
amplitude of the voltage on this capacitance. Eq. 
is rather general and is equally applicable to any type 
of conventional auto-oscillator (transistor, vacuum tube, 
tunnel diode, laser, etc.) in which the oscillation fre- 
quency is not strongly dependent on the auto-oscillation 
amplitude, i.e. in the limit du>/d\a\ 2 — > 0. 

There exist, however, auto-oscillators for which the 
oscillation frequency exhibit a strong nonlinearity N = 



dtjj/d\a\ 2 that is too large to be neglected. In such sys- 
tems, one expects that even small fluctuations in the 
amplitude at steady state can give important contri- 
butions to the phase noise. A pertinent example of 
present interest is the magnetic spin-torque (ST) nano- 
oscillator [1, & S H 0| ' which consists of a nano-sized 
metallic contact attached to a magnetic multilayer. Di- 
rect electrical current passing through the nano-contact 
can lead to a transfer of spin-angular momentum between 
magnetic layers in the stack [jjlll], which in turn creates 
an effective negative damping for the magnetization of 
the thinner ("free") magnetic layer. This negative damp- 
ing, analogous to the role played by an active element, 
can lead to self-sustained oscillations of magnetization in 
the free layer. The frequency of these auto-oscillations is 
determined by the applied magnetic field, static magneti- 
zation, etc., and is, in general, close to the ferromagnetic 
resonance (FMR) frequency in this layer. 

In contrast to traditional quasi-linear auto-oscillators, 
the frequency of the ST nano-oscillator strongly depends 
on the amplitude of the magnetization precession a: 
tola) — ujQ+N\a\ 2 . The sign and magnitude of the nonlin- 
ear frequency shift coefficient N depend on the direction 
and magnitude of the bias magnetic field (see 0, 0, 0] 
for details), and can be varied over range comparable 
to the oscillation frequency itself. Thus, the classical 
result ([1]) cannot describe quantitatively the generation 
linewidth in ST oscillators, and a new theory that ex- 
plicitly takes into account the nonlinear frequency shift 
of the auto-oscillator is necessary. 

In this Letter, we develop a theory to describe the 
generation linewidth in an auto-oscillator with nonlinear 
frequency shifts, and show that this nonlinearity loads 
to significant linewidth broadening. The theory is then 
applied to the ST nano-oscillator and we demonstrate 
that the correct treatment of such nonlincarities is essen- 



2 



tial for even the qualitative description of the nonlinear 
auto-oscillator. 

It has been shown previously 

H ES EH that the dy- 
namics of the dimensionless complex amplitude a (nor- 
malized by \a\ 2 = (M - M Z )/2M , where M is the 
length of the magnetization vector in the free magnetic 
layer, and M z is the projection of this vector on the equi- 
librium magnetization direction z; see [ij for details), 
which describes magnetization precession in an ST os- 
cillator, is given by 



■£ + i(uj + N\a\ 2 )a + T (l + Q \a\ 2 )a 
at 



(2) 



-r s (i - Q, 



\a\ 2 )a 



fn(t) ■ 



ujq is the linear oscillation frequency and N is the non- 
linear frequency shift coefficient (see Eq. (38) in [§]). T 
is the linear part of the positive damping and character- 
izes the equilibrium linewidth in the passive regime, and 
Qo is a phenomenological coefficient characterizing the 
nonlinearity of positive damping (see [12] for details). T s 
is an effective negative damping introduced into the os- 
cillating system by the active element (in the case of a 
spin-torque oscillator T s = al, where / is the bias current 
and a is the spin-polarization efficiency defined in Eq. (2) 
of (HI), and the coefficient Q s describes the nonlinearity 
of the negative damping (Q s = 1 for ST oscillator [§]). 
f n (t) is a stochastic term that accounts for the influence 
of the thermal fluctuations (noise) . The function /„ (t) is 
a white Gaussian noise with zero mean and second-order 
correlator 



(fn(t)f:(t r )) =2T P n S(t-t') 



(3) 



Here P n is the oscillator power at thermal equilibrium, 
i.e. (|a| 2 )r s =o = Pi- Note that the oscillator energy 
can be written as E(a) = /3|a| 2 , where the constant [3 
depends on the normalization of the oscillator amplitude 
a. Then, the equilibrium noise power P n can be written 
as P n = k^T/ (3, where fce is the Boltzmann constant and 
T is the absolute temperature of the system. 

We would like to stress that while Eq. @ is obtained 
for the case of a spin-torque oscillator, it can adequately 
describe an auto-oscillator of any nature under the influ- 
ence of white noise f n (t), provided that this oscillator has 
a nonlinear frequency dependence (ujQ + N\a\ 2 ), nonlinear 
natural positive damping To(l + Qo| a | 2 )j & n d nonlinear 
negative damping r s (l — Q s |a| 2 ). 

The stationary solution of Eq. ([2]) in the absence of 
noise (f n (t) = 0) can be easily obtained, 



where 



Po 



a(t) = y/P e 



c-i 



— iu)t-\-i<l 



QsC + Qo 



uj + NP = oj + N 



c-i 

Q S C, + Qq 



(4) 

(5) 
(6) 



(f> is a constant oscillation phase and C = Ts/ro is the 
supercriticality parameter. The solution Eq. (|4|) exists 
only for ( > 1. 

Sufficiently far above the threshold (i.e. for P ^> P n ) 
the solution of Eq. (0) with the noise term included will 
be similar to the noise-free solution Eq. in the sense 
that the oscillation amplitude will be close to the mean 
value of V7b, i.e. |a(t)| = y/Po + SA(t), \SA{t)\ 2 < P , 
and the phase </> will be a slow function of time. Substi- 
tuting the expression 



(7) 



for a(t) in Eq. @, and retaining only the terms of the 
first order in 8A, we find equations for the amplitude 
fluctuations, 



dSA 



2T eS P 5A = Re(f n (t)e- i *), 



(8) 



and phase fluctuations, 



dt 



2N^P~ 5A = -4= Im(/ n (i)e-^) . (9) 
v Po 



Here the effective damping T e g is expressed as 

r O ff = r (Q s C + Qo), 



(10) 



and f n (t) = f n {t)e luJt . Note, that the statistical proper- 
ties of fnif) and f n (t) are identical. Therefore, tilde will 
be omitted in the following text for simplicity. 

There is a significant qualitative difference between 
the behavior of the amplitude fluctuations described 
by Eq. ([8]) and the behavior of the phase fluctuations 
described by Eq. ([9]). Since the oscillation amplitude 
at steady state remains practically constant on aver- 
age in the presence of thermal fluctuations, \a\ w \fPo~, 
the correlation function for the amplitude fluctuations 
Ka{t) = (|a(i)||a(f + r)|) remains finite even if t — > oo, 
i.e. Ka — > Pa- Therefore, for large r the behavior of the 
full correlation function K(t) = (a(t)a*(t + r)) will be 
determined solely by the phase fluctuations, 



K(t) w P o (e l[0(t) - 0( 



t + T)}\ e 'iuJT 



fill 



For the frequency linewidth of the auto-oscillation, we 
are interested only in the fluctuations taking place in- 
side a narrow frequency region Auj -C r o ffPoi in which 
\dSA/dt\ - Au\SA\ < 2r eff P |<Jyl|. As such, the first 
(derivative) term in the left hand side of Eq. © can be 
neglected compared to the second term, and an explicit 
expression for 6A(t) can be obtained, 



SA = 



1 



2r off p 



Re(/„e~^) . 



(12) 



Substituting this expression for 5A(t) in Eq. §§§ leads to 
a closed-form equation for the phase fluctuations <f>{t) in 
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the system, 



at 




Re(/„e-^) + lm(f n e-^) 



/T5V 1+ fe |Im(/ " e »' (13) 

where a = atan(iV/r e gf). 

Eq. (fl~3|) is formally identical to the equation for phase 
fluctuations in a system without a nonlinear frequency 
shift (see e.g. second equation (9.8) in [l[), but with the 
increased noise level 



/»(*) - /;w = 




7n(*) • (14) 



Applying the general methodology to compute oscillator 
linewidths in the absence of nonlinear frequency shifts 
(see, e.g., Chapter 9 in [H or [3]) leads to an expression 
for the Lorentzian linewidth of the auto-oscillator with a 
nonlinear frequency shift N, 



Auj = r 



E 



N 

r ff 



(15) 



We have rewritten the ratio P n /Po as k&T/Eo, where 
Eq = (E(a)) = /3Po is the average oscillator energy. A 
comparison of the classical result (JTJ) with the general- 
ized Eq. (|15p shows clearly that the nonlinear frequency 
shift in the auto-oscillator leads to a significant linewidth 
broadening that is due to an effective renormalization of 
the phase noise (|14p by the nonlinearity N. 

The result in (fT5|) is the principal result of this pa- 
per and illustrates the fact that three key parameters 
determine the linewidth of an auto-oscillator with a non- 
linear frequency shift. First, the equilibrium relaxation 
rate of the oscillator To determines the overall scale of 
the possible linewidth variations. Second, the generation 
linewidth is proportional to the ratio of the noise energy 
(which increases with temperature) to the average energy 
of the auto-oscillation. Third, the ratio of the nonlinear 
frequency shift coefficient TV to the effective nonlinear 
damping T c ff gives a measure of the phase noise renor- 
malization due to amplitude fluctuations. 

For the particular case of an ST nano-oscillator an at- 
tempt to calculate the generation linewidth was previ- 
ously undertaken by one of the authors in but the 
nonlinear frequency shift was neglected. The calculation 
resulted in an expression for the linewidth (see Eq. (28) 
in [ijj]) that can be cast in the classical form ([I}, where 
the constant (3 is given by (3 — (Mo/jJuoVes, where 7 
is the gyromagnetic ratio, ujq is the oscillation frequency, 
and V a ff is the effective volume of the magnetic material 
of the free layer involved in the auto-oscillation. 

Note that the effective volume V c s for a ST auto- 
oscillator based on a magnetic nanopillar is equal to the 



volume of the nanopillar free layer itself, V c s = V, which 
is substantially smaller than the effective volume for a 
similarly-sized magnetic nano-contact. In the latter the 



effective volume is defined by Eq. (4) in [llj and is several 
times larger than the physical volume of the free magnetic 
layer under the contact due to the propagation of spin 
waves radiated from the nano-contact. Thus, our theory 
finds a natural explanation for a well-known experimen- 
tal fact (see [H, 0, H[) that the auto-oscillation linewidths 
associated with magnetic nanopillars are, in general, sev- 
eral times broader that those in magnetic nano-contacts. 

Another important result that follows from Eq. (|15p . 
in the context of a ST auto-oscillator, is the prediction 
of a linewidth minimum that follows from a change in 
sign in the frequency shift (e.g., from "red" (N < 0) to 
"blue" (N > 0)) as the magnetization is tilted out of 
the film plane. Across this transition the nonlinear fre- 
quency shift coefficient N passes through zero (see, e.g., 
Fig. 8 in [9j) at which one recovers the smaller value of 
the phase linewidth. A change in magnetization angle 
can be achieved in practice by applying a large exter- 
nal magnetic field at different orientations out of the film 
plane. This linewidth minimum has been recently ob- 
served in experiment (see Fig. 6 in 14]]). 

To illustrate the qualitative behavior of the generation 
linewidth of a ST auto-oscillator with the variation of 
external control parameters, such as the magnitude of 
the bias current / (which is proportional to the effective 
negative damping T s = a I in Eq. ([2"])) and the magni- 
tude H e and direction e of the external bias magnetic 
field, we present below results of the linewidth obtained 
from Eq. (fl~5|) . We consider a circular nanopillar oscilla- 
tor with the following parameters of the free layer: satu- 
ration magnetization ^qMq = 1 T, Gilbert damping con- 
stant ctQ = 0.01, nonlinear damping coefficient Qo = 0, 
thickness L = 5 nm, and diameter D — 50 nm. For sim- 
plicity we assume that spin-polarized current excites spin 
wave mode with a frequency close to the FMR frequency 
in a continuous free layer (see Eq. (37) in [§]). In that 
case the natural linear damping Tq and the coefficient of 
the nonlinear frequency shift N are given by Eq. (31) and 
Eq. (38) of respectively. 

While linewidth narrowing is expected with the in- 
crease in bias current due to the increase in auto- 
oscillation amplitude (Fig. QJi), the linewidth behavior 
as a function of the bias magnetic field is more complex 
than what might be expected from (|15[) . Indeed, one 
observes a non-monotonous variation in the linewidth 
around field magnitude and orientations that lead to a 
change in sign in AT. In such regions, there is a rapid 
linewidth narrowing followed by a broadening as the 
field is increased through the critical region in which 
the linewidth minimum occurs (9 e = 80° in Fig. [TJa). 
This complex linewidth behavior further is exemplified 
in Fig. [5] in which we present the role of the applied field 
orientation. The position of the linewidth minimum is 
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FIG. 1: Generation linewidth of a spin-torque auto-oscillator 
calculated from Eq. (|15|l as a function of: (a) negative current- 
induced damping T s = crl for fioH e = 1.2 T; (b) applied bias 
magnetic field H e for al — 1.0 GHz; for several angles 9 e 
between the bias field and the plane of the free layer. 




Field angle 6 (degree) 

FIG. 2: Generation linewidth as a function of applied field 
angle 8 e for: (a) three applied fields at constant a I = 1 GHz; 
(b) two bias currents at constant /j,oH e = 1.2 T. Inset of (a): 
Equilibrium linewidth To as a function of 9 e for noH e = 1.2 T. 
Inset of (b): Comparison of A/ = IS.u/2ti from Eq. (|15|l (solid 
line) , 10 x A/ from Eq. (dashed line) with 1 = 9 mA and 
HoH e = 0.9 T, and experimental data from Fig. 6a in Ref. Il4l 
(squares). 



clearly sensitive to the field magnitude and orientation, 
as it is the combination of these parameters that deter- 
mine the three critical parameters in (|15p . namely, To, 
r c fr, and N. This linewidth dependence on the field ori- 
entation is in good qualitative agreement with the exper- 
imental results from Ref. 14. A quantitative comparison 



of Fig. [2j3 in which the generation linewidth was calcu- 
lated for the experimental nanocontact geometry with 
the same parameters as above except for fi^Mo — 0.75 T 
and Q = 3. In contrast to the classical result (fl|), which 
predicts much narrower lines and a monotonous decrease 
in the linewidth as a function of 9 e , the renormalized 
phase-noise result (IT5|) gives a reasonable qualitative and 
quantitative description of the experimentally observed 
behavior, in particular, the linewidth minimum around 

e e re so . 

In summary, we have developed a theory of the gener- 
ation linewidth of a nonlinear auto-oscillator with a non- 
linear frequency shift which generalizes the classical re- 
sult (JTJ) . The additional nonlincarity in the oscillator fre- 
quency leads to a renormalization of the phase noise far 
above threshold. Applied to the particular case of a spin- 
torque nano-oscillator, the theory accounts for a number 
of characteristic, but previously unexplained, features ob- 
served in experiment: (i) General linewidth narrowing 
with increases in the bias current and oscillation ampli- 
tude (see Fig. 4 in Ref. [3); (ii) Presence of a linewidth 
minimum as a function of the external magnetic field 
orientation (see Fig. 6 in Ref. [13); (hi) Narrower lines in 
nanocontacts compared to nanopillars. 
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between the theory and experiment is shown in the inset 



